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♦ From Serat (1) 



SBiBOLS 


A ooemclent of ^ in differential equation of 

rolling. 

a .... height of wave, crest to hollow, inches. 

B .... coefficient of^^^in differential equation of 
rolling. 

B .... heaa of model. Inches. 

h .... block coefficient of model. 

M . transverse metacentrie radius. Inches. 

C .... a scale factor in the dlffei^ntial analyzer. 

D .... a scale factor In the differential analyzer." 

$t .... a scale factor in the differential analyzer. 

6 .... a scale factor in the differential analyzer. 

g .... acceleration of gravity, feet per second^. 

QM .... tx%nsverse metacentric hei^t, inches. 

az .... rioting arm, inches. 

.... coefficient of 0 in I’roude*s equation for ex¬ 
tinction of roll. 

Eg .... coeffiolent of in Froude’s equation for ex¬ 
tinction of roll. 

k .... radius of gyration of mass of model about a 

longitudinal axis through the center of gravity, 
inches. 

L .... length of model, inches. 

m .... midship section coefficient. 
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gear ratios in the differential 
analyzer. 

.... wetted surface, square inches. 

.... complete rolling period of model, 
starboard to starboax^, seconds, 

.... complete period of wave, seconds* 

•... time, seconds. 

.... phase angle; angle Toiy which inclina¬ 
tion lags wave slope. 

.... coefficient of j in rearranged 

differential equation used in integra- , 
tion. 

.... coefficient in rearranged dif¬ 

ferential equation used in integration. 
.... coefficient of QZ in rearranged diffeav 
ential equation used in Integration. 

.... displacement of model, pounds. 

.... transverse inclination of model, measured 
from vertical, radians. 

.... wave length, inches. 

.... height of raetaoenter above waterplane. 
.... maximum wave slope. 

.... slope of wave to horizontal. 

.... circular frequency of sine wave. 
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SUfiMARX 


A. OBJECf . 

The object of this investigation was to detei?- 
mlne whether or not the differential analyzer affords a 
practical means for predicting the rolling characteristics 
of ship models. 

B. MBTHOB. 

From data on a series of models tested by M.E. 

Serat (l), coefficients for the differential equation of 
rolling were computed, scale factors for the differential 
analyzer were calculated, and actual integration of the 
rolling equation was carried out with the differential 
analyzer, using a basic machine diagram as shown in Figure I. 
Results from the analyzer were then compared with the actual 
data obtained by Serat in tabular form (see Results, Table I). 


G. RESULTS . 

Actual plots of angle of inclinetlon from the vei>- 
tlcal as a function of time were obtained directly from the 
analyzer (see "Results", Figure II and Table II). 

D. OOHCLUSIOHS AND RECOMMEHMTIOHS . 

1. The method developed In this investigation of 
predicting the rolling characteristics of a model is funda¬ 
mentally sound. 

2. The method is also practical viewed from the 
standpoint of time required. A resume of time expended may 
be found in Results, Fart B. 
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In addition) the method is also practioal when 
viewed from the cost espeot, A resume of estimated charge 
for a fiin^e integration may he found In "Discussion of Re* 
suits" y ?art C. 



ItJTHODUOTIOIi 


She object of this investigation was to deter- 
aine whether the differential analyzer affords a prsotloal 
means for predicting the rolling characteristics of ship 
models. William Proude (2) proposed a method of solution 
of the differential equation of rolling by means of graph¬ 
ical integration. This method is quite lengthy and tedious 
suad as a result is rarely used. The differential analyzer 
installed at the Hassachusetts Institute of Technology was 
designed for the solution of differential equations and is 
readily suited to the solution of the rolling equation. By 
1320 methods shown herein, a curve of angle of inclination 
versus time can be obtained with a limited amount of given 
data. 

The Investigation was undertaken on a model only, 
because it was found impossible to obtain full scale roll¬ 
ing data after several attempts in different sources. 
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PROCSDUKS 

The essential steps followed in this investiga¬ 
tion were as follows: 

1, Evaluation of coefficients in the differential 
equation of rolling with the aid of a paper prepared by M* E. 
Serat (l). 

2. Evaluation of scale factors for the differ¬ 
ential analyser. 

3, Setting up analyzer and actual integration. 

4. Evaluation of results and coiaparlson with test 
data collected by Serat. 

For a detailed discussion of the procedure see 
Appendix A. 





y&Me X beXov shows a ooaparlson of Serat's ol>> 
aorreft data with that oalculated luidor the following eon- 
ditlona; 

San !• Som^-using ooiqplete rolling equation 
end SATe elope at the surface, 

San 2, i^dueed elope-uelng oompletm rolling 
equation and ware slope at a depth eorreponding to the 
dep^ of the oenter of huoyanof in still water. 

Ba& 3. Approxlaate-uslng rolling equation nodi- 
fled hf negleoting the danplng te» In § . 

Ihui h. 6onputed-from fozmtiae given In hossell 
and Chapsan 

Atm shown Is the "phase angle* representing the 
Wdiie hf which the an^e of Inclination lags the wave slope. 

yahle SI Is a representative seotlon of the 
recording tape of the differential analfser. 

yahle in contains the values of 9 versus tine. 
Figure XI Is a tracing of the plot nade bj ^e 
differential analfser during the ease ofole in idiloh the 
values in Shble SIX were recorded. 
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B* nm mmmB. 


Mn eatlBiate of time expended in obtaining the 
hljonre x^ealts follows: 

1. Svaluetion of coefficients 

in the differential eqnatlon of roll¬ 
ing 3 hours 

2. Staluation of scale factors 


and gear ratios in the differential 


analyser 

% 

flotting curves for input 

6 hours 

tallies 


3 hour s 

K 

Betting up analyser 

16 hours 

5 . 

Actual integration 

30 hours 

6. 

Evaluation of results 

H- hours 


Total tiaie eagE>end8d 

6S hours 


^ ABTmmmm mat , 

fhe coat involved in attaining the results s^s 
approxisatelf as follows: 

!• aeservation ehaxgee on dif- 
terentiaX analfxer, 3 ^ays at | 20.00 
per day t^.OO 

2 . Sourly chargee on integrators, 
hours on 3 units at $0*325 per unit 

per hour 22.5® 

3 . l*abor (bargee, hours for 

two operators at $1.60 per hour 76.S0 

#159-30 
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3PISCIIS3XOH OF mm*TQ 

A, yssaiOAi. . 

X« fhe aoQoracjT of results le such tl»t they oan 
bo aaauiftaA reliable to tenths of degrees so that there is a 
dlfferetioe of a half degree betsreen experiaentAI data and 
ilNai 1* ferat estiaated that the roll votxld be 16*0 degrees 
a aethod giren in his artiole. fhtzs the anal/ser gaTs a 
^^ser estiaate. Ihrors renainifig between 13*2 and l4,7 
degrees, or 3*^^ error, are probably due to errors in dai^ 
ing itself, inability to read damping ourre to that 

aoottraey, or a ooabination of both. !Biat part due to errors 

* 

in the analyser was United to 0.4^. 

2 . use of redueed eare elope yields the same 

error but in the opposite direotion so that no conclusiTe 
statement oan be aiade* Sowetrer, the method used for deter¬ 
mining the restoring force, (l.e. or wedge angled 

does not Justify using a slope less than that at the surfaee. 

3. fhe use Of a modified differential equation 

also yields nothing on vhloh to base a definite oonoluslon 

other than with damping terms of relatlye values as used it 

• 2 

is unnecessary to include the term in 0 • For other values 
this term probably oan become important. It is of interest 
to note that the phase angle was the same in Bans #1 and i3. 

4 . the phase angle relationship did not oheek at 
all well, and this is readily explained by the fact that the 


0 



l!OztiQl.a glTen In Jloeadll and Chapman ^3^ where x Is the 

~n 

j^i&ae angle *. x s _in- an 

. 7 - , -(^y 

ean be transformed to x = 

~r, a.GM-^ ^ 

S Ti^ 

Uhls formnla was derived fro® the assumption that rioting 
arm la a linear fnnctlon of O. In the oase analysed the 
rioting ar® at points above € degrees was considerably 
above the line whose tangent is as anoh as 6*7$ above 
at 16 degrees* Since the teras in the denominator of ( 2 ) 
are of the same order of sagnitnde, a change in dl of less 
than 6 ^ can ehsnge the tangent of the phase angle enough to 
give an x of the value attained above, or about 63 degrees* 

fhe record of time expended, given in ^llseults B*, 
does not, of eourae, refleot the true value of the method. 
1Sh0 authors at the outset were unfamiliar with the analyser 
and the aetlu^ds used in its operation. & person experlenoed 
in the use of the analyser would obviously require much less 
tins* Below is an estimate of the tine that ^ould be re¬ 
quired for a methodical attach on this problem, providing 
borough data has been made available: 

1 . Svaluetlon of coefficients 1 hour 

2. gvsluation Of scale footers 

snd gear ratios 5 hours 

3 . tlotting curves for inputs 3 hours 

4. Setting up analyser S houre 

5 . Sotual integration % hour per iwnrs 

6* Bvalimtion of results 1 hour 

fbtal time re<|uired 1 % hours 
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%t tills type ef probX» isere alimys vorlted on the analyzer 
Hie set-*t^ time would be reduoed fx^m S honre to about 3 
hourst and the eoale factor coopatations and plotting could 
be done in about 2 hours apiece* ^us» Hie total would be 
reduced to about % hours* 

gOSf * 

As wlHi tiae^ the cost is in excess of that which 
would hold for an lirrestlgator sore experienced with the 
probleB and this type of solution. An estioate of cost for 
an Integration for one ship or BOdel« and for three or four 
sea conditions^ is as followst 

X. Reservation oharge 120*00 

2. i^bor 25* 60 

3 . Integrator charge* 6.00 


fotal cost 


151.60 
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msamsiQMs 


A. 

fhe aethod beraln described does afford a prae- 
tloal aeafis of solving the differential equation of rolling 
both IToa the standpoint of tlae and cost, assoaing that 
the results obtained are considered of iaportanoe to the 
deeig^er. It would se^ that if oharaoteristies of this 
natuire were desired at ^1, this would be the logical ap-> 
proaidi. It is obviously sore direct, accurate, and less 
tedious than the graphical integration proposed by froude 

m. 


m 

It is also concluded tlmt for oases shere %2 
is less than of solution can be obtained with¬ 

in acouracy with the tem in 4^ onitted. It appears 
to have no effect on the magnitude of roll, velocity, ao- 
oelei^tlon, or even phase angle, fhon restricted to the 
value above. 



BEXKSMMEHD&TIOHS 


i? 


let l 0 felt that thie aethofi deserree further 
liiTestigatlQQ with the poeslhllltT^ of praotlcal uee In 
oonneetion with prellainarjr design. It does not seem Im-* 
praotioal that this type of solution could he accomplished 
on a special type of differential analyzer built for the 
specific purpose of solving only the rolling equation so 
that a naval architect could have his own apparatus. This, 
of course* would depend on the iaportanoe that the naval 
ardbdteot attaches to this more elaborate method of in- 
vestigating a design’s sea-keeping ability. JPerhaps some 
is»rt of standai:^ waves could be established to give a basis 
of comparison among several ships* much as is now done in 
longitudinal etrength. 

W, 

18ie extension of this method is hampered by lack 
of rolling data* and not a great deal of information of 
ship’s damping is available, fhese have probably been ne- 
^^ted heretofore due to the fact that there has been no 
particular use for such data. Sowever, if this method were 
put into practice* such data could be obtained end investi- 
gatlone made to correlate it to a store advanced point than 
is now the case. 



APPEHSXX 



x$ 


APfSHDIX A 
mthliS OF gaOOEPgRS 

X» ^alnatico of coofflotentg i n the differential egnation 
of rolling . 

fhe differential equation of rolling used for in- 
teipration iraas 

+• A0+■ S©*‘+- O ^1) 

% 

Before integration it vas neoessary to evaluate g» i, 

% and ^e value of g was taken at 3^^ inches per (seo- 
She paper of Serat Cl) gi^ee information on several 
BOdele of the parallel middle bodies of ships. One of these 
models*, il -2 with bilge keels Installed, was selected for 
study. Its idaersoteristics are given in Figures III, I?, 
and f, and fable I?, the value of was given but other 
tezms were found by less direct aelhode. In order to get & 
and B it was necessary to eonvert the curves given by ierat 
of the rate of angular damping versus angle of roll for 
rolling in still water to a curve of extinction. This was 
dona by multiplying the rate of angular damping at a given 
angle by the period of roll in still water and plotting this 
decrement against the angle of roll, laving the curve of 
extlnotlon, its equation was approximated by: 

b© = K.©+-(2) 

and values of and evaluated. It so happened that 

the curve could have been approximated by the formula; 
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FI6UR6 HL 


GENERAL FORM OF MODEL n-2 




TABL-tlS: 

MODEL DATA 

DISPLACEMENT 10.21 LB. 

LENGTH 20.00 IN. 

BEAM 6.64 IN. 

DRAFT 2.25 IN. 

MIDSHIP-SECTION AREA 14.14 SQ.IN. 

MIDSHIP-SECTION COEFFICIENT 0.96 

WIDTH OF BILGE KEELS 3/16* 

WETTED SURFACE 201.6 SQ.IN. 

HEIGHT OF METACENTER ABOVE W.L. 0.67 IN. 

METACENTRIC RADIUS 1.76 IN. 

BEAM-DRAFT RATIO 2.91 



























uc.\-«rvc-rN^HN I pvwu-l- 


Cl 


FIGURE EC 


DAMPING CURVES MODEL 11-2 
WITH BILGE KEELS 
GM *0.496" T* 1.18 SEC. 
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eaa» in use of the differential analjrser stade it 
ffiore feaslhle to use the first equation girea. With the 
aid of the follotring equations from Sossell and Ghspman l 3 l*» 
the Talues of A and B were obtained once and Eg were 
found. 

K 

I3) 

v<r = 

. C45 

In these equations the period of resisted rolling in still 
eater, 1 ^, and setacentrio height, ears given Serat trm 
hie still eater rolling test. 

Also froa Ej, fj,, and elth foraulae given in Bossell 
end Ghapaan 1|3J* value of k was found 

».ioS-<fe I 



Einalljr, the righting aiia, S2, la not a constant in 
the Integration and wae considered to he a function of the 
iaaersed wedge angle shioh is the difference between the angle 
of inolinatiott from the vertical, and the nornal to the 
wave slope, ff* as shown below: 




2k 


the etatieaX stability mvwe^ versus #, in 
etiXl water which is aotuaXljr rightii^ arm vermis wedge 
angXe^ was need as an Input to the analyser with the arga^ 
■ent 9 - f, and referred to hereafter as t(9 - 

Wave elope, fi, was ooneidered a sine funotion of 
tiae and the maxlmm value «ae eoaputed frost 

radians Ibl 

a curve of p versus t vas used as an input to the 

analyzer. 

Svaluation of scale factors and ^ear ratios . 

Since the differential analyzer ie a aeohanical 
apparatue It has definite physical limits. In order not to 
overrun these limits, and in order to use them to beet ad«» 
vantage, estimated extreme values of ell variables involved 
must be eeleoted and scale factors assigned to all variables. 

A diagram such as figure TI following was di^wn up 
and extreme values estimated and called 

fmf %i» K* 

fhe differential equation {it mu tx^sposed so that 

0 + c<6 4-^©^ + ■^-P(0-cp)-= O 

viiere 

fteiurii^ and scale factors were denoted ty letters 
me i^wn, and their talues were eosimted fzom oertaiit ex¬ 
pressions that can be written, flnraughout the oaloulatlons. 


o< = :5s. • >r- ^ 






aesld t&otora rater to the nxmioer of tums of a gi^tn shaft 
to pspreseat on© unit of the variable. For exeaple, the 
tiae dteft ©akes A turns to represant on© unit of ti®©* 

0ear ratios ar© given as the ratio of thriven shaft turns to 
driven turns. Integmtors produce the product of the soals 
factors of the Inputs divided W 32. For example, Integt^tor 
#1 takes mfii and ^t, and produeas VTith these facts 

32 , 

"knowa^ all veriables 'with appropriate scale factors were 
written at the end of each ahaft as shown. Ih© ©qualities 
and inequalities were set up as followss 

Ca) Since the pOler plot of p versus t makes on© 
turn for on© hundred turns of the input shaft and the period 
of the wave was known 

HjAt y^Ti =■ 100 

(b) Since the polar plot output follower he© a 
Hsiit of ^ JJIO turns fro® midposition 

8-cp^<'80 

Cc) Input plot 2 -me on a table with a rang© of 
the independent variable of ^1^0 turns, and of the dependent 
variable of 36 O turns, ilenc© 

n^5'<2(9-(pL«4.80 
D* 2[V+(e-(p^J <ZbO 

(fhe factor 2 resulting from the fact that both variables 
sKiy be plus or minus,) 


< 10 ) 

(11) 
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(d) the integrators have an excursion of 40 turns 
to either side of their raidposition, hence 


n^ce„ < 40 

(W 

ntoD©« < 40 

(1J> 


te) ^e equalities fro® the adder gears ares 


'OtiVyg AC 

im 

<oA- ^ ^ 


A.T) • 

( 16 ^ 


Solving these equations by trial and error, the 
gear ratios as shoen in figures fll were arrived at. 

Searing for record of © versus to and their 
counters ms arranged to give valueo easily read, as 10©, 


looot, etc. 

3 * Betting up analyser and actual integration . 

ftie analyzer ms arranged as shown in figure Til, 
and the first integration, hereafter referred to as tun 1, 
was carried out with the wave elope ©f the surface, fhe 
maxiamm angle of roll in the steady state and the phase re¬ 
lationship between wave slope, or emiting force and result¬ 
ing oseillatioa ms estioated frost 



whiOh are given in Hossell and Chapoan C3K 
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" SCHEMATIC DIAGRAM 

differential ANALYZER 5ET-U? 

UE 6 END 
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4. PROI^ UV5H GEAR — 1 >- 7 . RjOt-R -3— 

ARROWS INDICATE FOSITIVE: ROTATION A5 VieWEi? FROM ABOve- 
FIGUR6-S ON SHA1=T tNDlCAre NUMBER OFTURN5 FDR I UNIT OP VA«a-l4«k&u& 
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With these Talues initial oondltlone were set np 
vhere t vas sero» angle of Inollnatlon was # eoc^uted froai 
Cl7)» the Teloolty of the roll was aero, idle ware elope was 
leading the angle of roll by x degrees and the aocelez^tlon 
was the n^atlve of f(0 *- fhe problem was run until 
differences between successive inelinations were less than 
for several conQilete cycles, fhie was assumed to be 
the steady state. Tapes recording counter readings approxi¬ 
mately 50 tinea per cycle and record of © versus t plotted 
by the ansdyzer were made during the run and were marked when 
the steady state was reached. 

Following the same general procedure Run 2 isisis 
made using a wave slope corresponding to motion of wave par¬ 
ticles at the height of the center of buoyanoy. This was 
done by oonsiderlng that the wave id.ope decreases as an ex¬ 
ponential fi^m the surface so that 

-Wn-ct , 

Cp'- ^ -TX (19) 

where d is the depth below the surface to the center of buoy¬ 
ancy and ^ is the maximum value of the wave slope at the 
sarfaoe given in equation (6). 

The final integration^ Run 3 , was made with the 
sarfaoe wave slope plot, but with the #2 Integrator inopera¬ 
tive, so that 

0-4-o<©+ 'ZJ'f* "’tP) “ 0 

was being integrated. This was done to show the effect of 
the © term, which is often negleot^. 
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Xfl ftXl runs the i^aee eisgle at steadi’ state mB 
iTedordled at the end of the run. 

i^aluation df results . 

I^e plots of 0 versus t were used oaljr as a guide 
to the general trend of the Integration and results were taken 
from printed record of counter readings. Counter readings or 
turns nade were conrerted to the actual values of the vari¬ 
ables by application of the appropriate scale factor. Maxl- 
SRUI values of inollnation were compared with the estimates 
and etperiaental values given by Semt (iK 

A record was prepared of time spent on various 
phases of the work so as to contrast this with the other 
Betlu»ds of estimating or oomputing roll to be eapeoted fzom 
a ship fomt among waves. 
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APPENDIX B 


SAMPLE GALCULAflOHS 


Evaluation of coo^fictents * 

The curve of angular damping in Figure IV was con¬ 
verted to a curve of extinction multiplying I.IB seconds. 


Nhen 


e iz* 


50= S.S'zVroW 
50= 8.'4-Z"/^oH 


50= 0.^2.V.r« 

By equation (2> 

l2K,+ »44-Kt= -S.S2 
Kbi<,'«-2S(b S.4-2 

Solving simultaneouely 

K,= 0.2(b2 l<2,= O.OlfoS 

iFroa equations (j) and ^4) 


C»o-2i)(49(al_ 0.»S8 \k)-i*^-:5ec 


(.0>faSVl0.2» V.< 


_ 0.0022 I to- iV> - -,ec^ 


equation (5J 


ife --LLSdJbik. 


_ (.2fa-2)^ - 2.53 IKI 


o< = Ag _ (.>S6)(38(b) _ 0.693 /4ec 

(I o.- 2 1)02.69)^ 


/3 = = (.ooeaYsSfa) ^ o.o»2A 

A-fe*' C'0-2»)(2.59Y 


« 3SCg_ ^ 
" (2.S9)^ 




/n oe:: 
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Bf eifuation and ixaing irare of 1.3f inohes 
hei^t mA 9^ Inches Xength 

U)^ = o. 04&4- »«!«»» 

' or 2.67® 

Wov the m&te of reduced slo|>es 
fron t&hle Zf A « 1.1^ Ijsehe# 

From equation <X9l 
-z-jrA-ito). 

(p^ _ e ■ S^(.o4-<o4)= o. 04^*0 radians 

or 3.47® 


Evaluation of sogJle faetors and gear ratios . 
Baximugi value of # eaiieoted: 

Q^^±2.0^ or ±0.2»Sr-cxoL. 




Jf #^^c5«ra>f irlth period of 1.^ seconds 

I ©<vn I ~ ©AV» 

I ^/vhI ~ 


ccnoC 00 - -§^ = 5.24- 

therefore , • , 

l.ft3 rotOL/sec 
S.SS" rrxcLI&zc^ 

'*** Gm.- -O.'SSra.ci. ‘fU= -.O^DrocoC 

=3 - 2.S rxxdy^bec 

0/m = - t3.ara.ci/5)ec* 

At (0 ts t .41, from stability ourte. 

= ±0«3€>in 





Sabstitutlng the so in ecjuatlono W throng {10 
and hgr txlal and error the following were established: 


As 1^29 
» s 2^90 

e s 152 .% 

M s 


n.- 0 04S5^5---^(H^i) 

/fz - 
fl^ =r 

ttjf. 3 / 32 . 

dir- 'A 
d6 - V3 
d; " YiozA 


Getting OP analyzer and actual integration, 

Sinee the ware slope ims considered to be a siO'- 
usoidal function of time the eurre on titie polar plot was a 
oirele, this plot rotated once for each l*2(ld29) turoa of 
the t ihaft^ thus representing a sine function with a period 
of 1,20 seconds idnee 1S29 turns were equiwalent to one sec¬ 
ond, l^e output was 24^0 fo represent a maxlmuai wave 
slope of 0,€^C4' radians would take 2450(0,0464J turos. 

With a follower making 20 turns per inch the circle should 
be made to hawe a radius of 

^ss^mia ^.sr tn^B 

the input on table #2 was plotted by the follower 
to be used in the runs in order to awold errors Inherent in 
the lead screws. Points were selected from Figure W, oon- 
werted to turns by 612,6(0 - p) and i K tip - PK 




Sfttlaated ttaximm angle of roll and phase angle 
were frea equations (171 (l®)- 


e = 


^ 0.O4fe4- - ^ 0.263 


l.-zo*- 




X = -/■o.K.-' —L^o_~ '■hx^-‘^e> 

^ '-i¥ioT 

x= 76®20' 



APPSHDIX 0 


OaiaiHAL 

I^ta as aetoalljr recorded bjr the analjrzer is 
filed with the Center of Analxsis, Sassaohueetts Institute 
of fedhnology. 
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